Mapping the phase diagram of strongly interacting matter 
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We employ a conformal mapping to explore the thermodynamics of strongly interacting matter 
at finite values of the baryon chemical potential fx. This method allows us to identify the singularity 
corresponding to the critical point of a second-order phase transition at finite given informa- 
tion only at = 0. The scheme is potentially useful for computing thermodynamic properties of 
strongly interacting hot and dense matter in lattice gauge theory. The technique is illustrated by 
an application to a chiral effective model. 
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INTRODUCTION 

The properties of strongly interacting matter at finite 
temperature, T 7^ 0, and baryon density, ub 7^ 0, have 
been studied extensively in recent years. The phase dia- 
gram has been explored experimentally in heavy-ion col- 
lisions and examined theoretically in models as well as in 
first-principle calculations of Quantum Chromodynamics 
(QCD). In the nonperturbative regime of QCD, Lattice 
Gauge Theory (LGT) provides a powerful tool for com- 
puting the thermodynamic properties of strongly inter- 
acting matter. 

In LGT calculations, the Feynman path integral for 
the partition function of QCD in Euclidean space-time 
is evaluated using Monte-Carlo sampling. The quark de- 
grees of freedom are integrated out, yielding an effec- 
tive action for the gluons, which includes the so-called 
fermionic determinant. For nonzero (real) values of the 
baryon chemical potential /i, the determinant is complex, 
which implies that the Monte-Carlo weight function is 
not positive definite and consequently that this technique 
fails for Re /i 7^ 0. This is the so-called sign problem, 
which has impeded progress in lattice calculations at fi- 
nite baryon density. 

Various indirect approaches have been developed, to 
sidestep the sign problem (for a review see Ref. [1]). 
These generally involve an extrapolation from ensembles, 
where the fermion determinant is positive definite, e.g. 
from /i = or from imaginary /i. Using such techniques, 
thermodynamic functions can be computed for small real 
/i/T, i.e., for small baryon densities. 

The extrapolation from the /i = ensemble can be 
performed, e.g. by means of a Taylor expansion [2] in 
/i/T. Without refinements, this method is applicable 
only within the radius of convergence of the series, R^. 
Also for imaginary /i, the fermion determinant is strictly 
positive definite and consequently systematic LGT cal- 
culations are possible An analytic continuation to 
real /i by means of a polynomial is valid only within the 
convergence radius of the Taylor expansion Thus, 
both methods are restricted in applicability by R^. 



The radius of convergence of the Taylor series is limited 
by the distance to the closest singularity in the complex /i 
plane. Conversely, the convergence properties of a power 
series yields information on the singularities of the orig- 
inal function. Of particular interest are those singulari- 
ties corresponding to the critical point of a second-order 
phase transition, to a crossover transition 0, [11 or to a 
spinodal line 0- We also consider the "thermal" singu- 
larities associated with zeros of the inverse Fermi-Dirac 
function [9]. 

The convergence properties of the Ta ylo r series in /i has 
been studied in model calculations [9|,llo|- It was found 
that on the order of 20 terms or more are needed to ob- 
tain reliable information on the structure of the phase 
diagram for this model. This is well beyond what is 
presently available in LGT calculations for small quark 
masses llll,[i2|. Here more refined methods for the anal- 
ysis of a truncated power series, based e.g. on the Fade 
conjecture (for an application to QCD in Ref. (HI) or on 
the method presented in this paper, may prove useful. 

In this article we present a method for unravelling a 
singularity in the complex /i plane, connected with a 
second-order phase transition, given a finite number of 
terms in a series expansion of some thermodynamic func- 
tion. The method, which utilizes a conformal mapping 
of the Taylor expansion in /i, yields reliable results al- 
ready at fairly low orders of the series. Such methods 
have been employed in quantum field theory for ana- 
lytic continuation of perturbative results to the strong 
coupling regime Q, as well as in condensed matter 
physics m, [TgI, where conformal mapping techniques 



were used to enhance the sensitivity to the properties 
of the critical point [l7j. 

Methods for analyzing a power series of the type em- 
ployed in this paper, are based on the theorem of Dar- 
boux on asymptotic approximations to the late coeffi- 
cients of a series [l8|. These methods can be used if the 
physical singularity, e.g. critical point of a second-order 

If this is not the case. 



phase transition, is dominant 
it may be possible to transform the series by means of a 
conformal mapping into a form where the critical point is 
the closest singularity, thus permitting an analysis based 



on Darboux's theorem. Such a transformation enhances 
the sensitivity of the series coefficients to the physical 
singularity and minimizes the influence of other singu- 
larities. 

In another class of methods, Pade approximants are 
used to improve the convergence properties of the series. 
These methods, which are based on the Pade conjecture, 
are, in principle, not restricted to the case where the 
physical singularity is dominant. Nevertheless, the re- 
sults of such an analysis are usually more reliable also in 
the case of one dominant singularity [l^- In this paper, 
we explore a method based on Darboux's theorem, but 
stress that methods of the Pade class also profit from a 
judiciously chosen conformal mapping, which enhances 
the sensitivity to the critical singularity. 

The analytic structure in the complex /i plane was also 
explored in a simple exactly solvable model in Ref. 0, 
where in particular the role of Lee- Yang zeros in finite 
systems was discussed. In this article, we focus on sys- 
tems in the thermodynamic limit, and defer the discus- 
sion of finite-size effects to a subsequent paper. 
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CONFORMAL MAPPING 

The thermodynamic functions are analytic in the cut 
complex /i plane. The cuts are associated with the sin- 
gularities mentioned in the introduction. First let us as- 
sume that the location and type of the singularities of 
a thermodynamic function in the complex /i plane are 
known, but only a finite number of terms in the Taylor 
expansion about /i = can be computed. In this case one 
can employ a suitable conformal mapping to improve the 
convergence properties of the power series and thus ob- 
tain reliable results beyond the radius of convergence of 
the original expansion [28| . 

To illustrate the idea we consider an elementary ex- 
ample, the chiral quark model [2l| (QM) in the mean- 
field approximation. The method is, however, general 
and applicable to a wide range of problems. In the chi- 
ral limit, for vanishing chemical potential, /i = 0, the 
system undergoes the second-order phase transition at 
the temperature T = T^. In the QM model, the criti- 
cal temperature decreases with increasing chemical po- 
tential. Then, at a temperature T < (and larger 
than the tricritical temperature), a critical point is lo- 
cated at some value of the chemical potential /ic(^) 
and by symmetry at —fic{T) (see Fig. [TJi). Owing to 
the periodicity of the Fermi-Dirac distribution function 
fp{T^jii) = l/{exp[(£p — lJi)/T] + 1}, the branch points 
are repeated in the imaginary ji direction with the pe- 
riod 2z7rT. 

In addition to these singularities, there are also thermal 
branch points, corresponding to the zeros of the inverse 
Fermi function that are closest to the imaginary axis [29| , 
/ith = ±[min(£p) +i7rT + 2z7r/cT], k eJ\f. In QCD, owing 



FIG. 1: a) The analytic structure of the order parameter in 
the complex /x plane. The branch points at di min(£p) zb mT 
are the thermal singularities (open dots), while those at ±/ic 
on the real axis correspond to the critical point of the assumed 
second-order phase transition (filled dots). 

b) The analytic structure of the order parameter in the com- 
plex fugacity A plane. The thermal cut is located between the 
branch points at and — oo, while the branch points at Ac and 
1/Ac are associated with the second-order phase transition. 

c) The analytic structure of the order parameter in the com- 
plex w plane. All cuts are on the circumference of the unit 
circle. The branch points corresponding to the critical point 
are at it; = 1 and = — 1, respectively. 

d) The analytic structure of the order parameter in the com- 
plex w plane for Ag > Ac. The branch points corresponding 
to the critical point are located w = Wc and w = —Wc, 
respectively. 



to the Z(3) symmetry (see Ref. [3]) the thermal singu- 
larities are located at Im /i = ±[7rT + 27r/cT]/3. 

A more transparent picture emerges when we apply the 
mapping A = e^/-^, where both the thermal cuts and the 
copies of the critical branch point are mapped onto the 
real axis. This transformation amounts to a change of 
variables from the chemical potential /i to the fugacity A. 
In the complex fugacity plane (see Fig. [TJd) the thermal 
branch points at juth are mapped onto the negative real 
A- axis, while images of the critical points are located at 
Ac = e^^/^ and at l/Xc The thermal singularity closest 
to A = 1 is the branch point at Ath = 0, corresponding to 
Ep = oo. This is a universal property, valid also beyond 
the mean-field approximation. Consequently, since Ac > 
1 and hence < 1/Ac < 1, the singularity closest to 
A = 1 is the one at l/Xc Thus, in the complex A plane. 
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the radius of convergence, = 1 — 1/Ac, is not affected 
by the thermal singularities. 

The critical points at Ac and at 1/ Ac are branch points 
in the complex A plane. In the mean- field approxima- 
tion, the order parameter scales as cr(/i) ~ (/ic — m)^^^ 
for ^ /ic- Thus, the analytic continuation of the or- 
der parameter exhibits one or more cuts starting at the 
branch point /i = /ic- In systems of finite volume, these 
cuts correspond to lines of singularities, the Lee- Yang 
zeroes 0,11, 22-111, which generally form a nonzero an- 
gle with the real axis. In the thermodynamic limit, we 
use the freedom to choose the location of the cuts and 
place them on the positive real axis starting at Ac and on 
the real axis in the negative direction starting at 1/Ac, 
respectively. 

We now map the cut fugacity plane onto the interior 
of the unit circle. A suitable conformal mapping is given 

by 



w{\) 



VAAc - 1 - yxT^ 

VAAc - 1 + VAc - A ' 



(1) 



The branch points at Ac and 1/Ac are mapped onto w = 1 
and w = —1 respectively, while the corresponding cuts 
are mapped onto the circumference of the unit circle, 
as illustrated in Fig. [Th. We note that the symmetry 
/i ^ — /i corresponds to w ^ —w in the mapped vari- 
able. A Taylor series about w = 0, which corresponds 
to A = 1, converges for all points within the unit cir- 
cle. Consequently, one can use this mapping to construct 
an analytic continuation of the Taylor series in A, which 
is vahd also beyond the radius of convergence in the A 
plane, R^. 

We consider first the ideal case, where the dependence 
of the order parameter on the chemical potential, or 
equivalently on the fugacity, is known. Then, using the 
inverse transformation to ([T]) 



A: 



Ac(l-^)^ + (l + ^)^ 

(l-^)2+Ac(l+^)2 ' 



(2) 



we obtain the order parameter as a function of 
(^wi'w) = cr{X{w)). The corresponding Taylor series in 
w is given by 



where 



1 d"a^ 
n! dw^ 



(3) 



(4) 



w=0 



The conformally mapped expansion as a function of A, 
which converges beyond i?^, is obtained by introducing 
the mapping ([T]) in the Taylor expansion (j3j) 



<Tcm(A) = <7»Mi)) + ^ <MA)r . 



(5) 



The convergence properties of the expansion ([5]) are im- 
proved compared to those of the original expansion in A 
or /i, because information on the analytical structure of 
the function is included in the mapping. 

RESULTS AND DISCUSSION 

In applications, e.g. to lattice QCD, a strategy by 
which one can find the location of singularities in the /i 
plane, given a finite number of terms in the Taylor expan- 
sion, is needed. We have found the procedure presented 
below useful for locating a second-order phase transition, 
i.e., a critical point on the real /i axis. In the case of QCD, 
this could be a critical endpoint, where the crossover 
transition ends and is replaced by a first-order transition. 

First we define a map of the type ([T]), replacing the 
critical fugacity Ac by a parameter Ag 



^g(A;Ag) 



yAA^-^A^ 



A 



(6) 



n=l 



\/AAg - 1 + ^ 

The analytical structure of the thermodynamic function 
in the Wg plane now depends on the value of Ag relative 
to Ac- For Ag > Ac, the branch points associated with the 
critical point, Ac and 1/Ac, are mapped onto points inside 
the unit circle at = iiL'c = ±tt;g(Ac;Ag) (see Fig{T]d). 
Since w = ±k;c are the singularities that are closest to the 
origin, the radius of convergence of the Taylor expansion 
in w is given by R^ = Wc. If, on the other hand, Ag < Ac, 
the critical point is mapped onto the circumference of the 
unit circle and consequently the radius of convergence 
equals unity, R^ = 1. Thus, given R^ as a function of 
Ag, the location of the critical point, Ac, can be obtained 
e.g. by applying the inverse mapping to w = R^ . Hence, 
the dependence of the radius of convergence on Ag can be 
used to determine the location of a second-order critical 
point. 

In reality the radius of convergence is known only ap- 
proximately, because only a finite number of terms of the 
Taylor expansion are known. In order to enhance the sen- 
sitivity to the location of the critical point and minimize 
the influence of other singularities, it is advantageous to 
use a mapping which leaves the singularity of interest 
close to the origin and moves all others as far away as 
possible [13] • As we show, this can be achieved rather 
efficiently by varying Ag in Eq. (|6]). Let us assume that 
the Taylor expansion in /i or equivalently in A is known 
up to n-th order. We proceed by performing the map- 
ping of the truncated series and obtain an expansion of 
the type ([3]), truncated to n-th order in w. We note that 
the coefficients cj^, are uniquely determined by terms of 
order m < n of the original expansion. 

The radius of convergence as a function of the param- 
eter Ag is approximately given by R!^ = |c^|*^~^/^^. (In 
the limit n ^ oo this expression for the radius of con- 
vergence is exact.) The approximate expression is, for a 
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FIG. 2: The radius of convergence in the w plane as a 
function of 1/Ag calculated at fixed order n — 12 (dots). The 
error bars show the deviation of the radius of convergence R12 
from that obtained with n — 10, Riq. The solid line show the 
result of fitting with the Ansatz i?^(Ag). 
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FIG. 3: The phase diagram of the QM model in the chiral 
limit. The solid line represents the critical line of the second- 
order phase transition, while the dashed and dash-dotted lines 
are obtained using different estimates of the radius of conver- 
gence based on the original Taylor expansion in /i. The sym- 
bols show the phase boundary obtained using the mapping 
technique discussed the text. 



given n, fitted with the Ansatz i?^(Ag) = a + '^g(Ac; Ag) 
for < 0.95. The parameter a accounts for corrections 
due the truncation at a finite n, while Ac is the approx- 
imate critical fugacity. We neglect the dependence of a 
on Ag and justify this approximation a posteriori. 

The resulting fit is illustrated in Figl2]at a temperature 
T = 150 MeV for n = 10 and 12. By following the proce- 
dure outlined above for temperatures in the range from 
100 MeV to Tc, we obtain an approximative phase dia- 
gram of the QM model (see Figl3]). For comparison, we 
also show the results obtained by extracting the radius 
of convergence directly from the original Taylor series in 
/i, for n = 8, 12. The present approach yields a better 
estimate for the location of the critical point for all val- 



ues of n considered. The fact that over a wide range 
of temperatures, we find approximately the same phase 
boundary for different n confirms the consistency of the 
Ansatz for R^ . Thus, at least in the QM model, the Ag 
dependence of a is indeed very weak. The parameter a 
vanishes in the large n limit as expected, with the leading 
term proportional to 1/n. 

We have tested this method on trial functions with 
various types of singularities. In every case it was found 
to provide an efficient estimate for the location of the 
singularity. This indicates that our method is applicable 
to a broad range of problems involving the analysis of 
truncated power series. 

Another interesting method for locating a critical 
point, based on Darboux's theorem for the asymptotic 
expansion of late Taylor coefficients (see e.g. 118|), was 
proposed by Hunter and Guerrieri [26^ . This method has 
the advantage that it can be used to identify both the 
critical point and the critical exponent. We find that for 
a fixed critical exponent, P = 1/2, results of similar qual- 
ity as those shown in Fig. [3] are obtained, while attempts 
to determine also the critical exponent were not success- 
ful for n < 12. A detailed study of the two methods will 
be reported in a subsequent paper. 

For a physical pion mass, we expect a crossover transi- 
tion at small values of the chemical potential. In this case 
the critical singularity splits into two conjugate branch 
points, located at complex values of the chemical poten- 
tial symmetrically with respect to the real axis 0, Hll , /ico 
and /i^Q. Therefore, for a crossover transition, the radius 
of convergence of the Taylor expansion is given by the dis- 
tance of fico to the origin, \fico\- In the QM model, the 
resulting radius of convergence decreases continuously, as 
the temperature is increased from the critical endpoint 
to the pseudocritical temperature at /i = and then in- 
creases again. The minimum of the radius of convergence 
is thus located close to the /i = pseudocritical temper- 
ature and is not associated with the critical end point, as 
suggested in Ref. 



27|. 



In this article, we presented an efficient scheme for lo- 
cating the critical point of a second-order phase transi- 
tion in the thermodynamical limit. Generalizations of 
our approach to finite systems and to crossover transi- 
tions are under study. 
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